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Abstract 

Northern hemisphere TeV gamma-ray observatories such as Milagro and Tibet AS7 have demon- 
strated the importance of ah-sky instruments by discovering previously unidentified sources that 
may be the PeVatrons producing cosmic rays up to the "knee" in the cosmic ray spectrum. We 
evaluate the potential of IceCube to identify similar sources in the southern sky by detailing an 
analytic approach to determine fluxes of muons from TeV gamma-ray showers. We apply this ap- 
proach to known gamma-ray sources such as supernova remnants. We find that, similar to Milagro, 
detection is possible in 10 years for point-like PeVatrons with fluxes stronger than several 10^^^ 
particles TeV~^ cm~^ s~^. 

PACS numbers: 
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I. INTRODUCTION AND SUMMARY OF RESULTS 



The IceCube neutrino detector is projected to observe 220 atmospheric muon neutrinos 
per day in a background of cosmic ray induced atmospheric muons traversing the detector 
at a rate of 1650 Hz. Upgoing muon tracks initiated in or near the detector by neutrinos 
that traverse the Earth must be separated from a background of down-going cosmic ray 
muons, which dominates the neutrino signal by a factor 10^ at TeV energy. IceCube has 
demonstrated this capabihty with a partially deployed detector and has thus become the 
largest neutrino detector observing the northern sky as well as the largest instrument de- 
tecting muons from the southern sky In light of the recent discovery of sources of TeV 
gamma-rays that may be the PeVatrons powering the cosmic rays up to the "knee" in the 
cosmic ray spectrum 2i], it is important to revisit the potential of IceCube as a gamma-ray 
detector similar to the Milagro and Tibet AS7 arrays but located in the Southern Hemi- 
sphere. This paper will focus on the possibility of identifying PeVatrons in the southern sky, 
as well as known bright gamma-ray sources such as the supernova remnant Vela Jr. 

Recently, the Milagro and Tibet AS7 TeV gamma ray experiments have demonstrated the 
importance of constructing all-sky telescopes complementary to the more sensitive pointing 
air Qerenkov gamma-ray telescopes, discovering a class of TeV sources that had not been pre- 
viously distinguished by the pointing telescopes 3|, The properties of these non-thermal 
sources are nevertheless striking, with spectra consistent with an energy dependence 
possibly up to at least 100 TeV. They are the first candidate sources to be associated with 
the PeVatrons that accelerate cosmic rays to the "knee" at 3 PeV. The sources cluster in the 
direction of the nearest spectral arms, consistent with the idea that cosmic rays are acceler- 
ated by the remnants of supernovae exploding in star-forming regions. If these sources are 
PeVatrons, the TeV gamma-rays are the decay products of neutral pions produced in the 
interaction of cosmic rays from the remnant with atoms in the interstellar medium near the 
source. Supernovae associated with molecular clouds are a common feature of associations 
of thousands of OB stars that exist throughout the Galactic plane and so the flux detected 
at Earth should be larger from the vicinity of known molecular clouds where the greater 
target density is the source of an enhanced TeV flux. Observations are consistent with this 
scenario: Milagro observes gamma-ray emission with an average energy of 10 TeV over a 
broad patch of the sky in the direction of the star forming region in Cygnus as a background 
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to individual sources, which could be associated either with molecular clouds or with the 
supernova remnants themselves. 

We discuss here how it may be possible to identify similar sources in the Southern Hemi- 
sphere by operating IceCube as a gamma-ray telescope. The possibility of using muon 
detectors in general, and IceCube in particular, as a gamma-ray telescope has been exten- 
sively =tudwflflBay.yfl »jd has been exploited to put IMts on the TeV enh^ion 
of a soft gamma-ray- repeater [12|, Il3| . Showers initiated in the atmosphere by TeV gamma 
rays from sources in the Southern Hemisphere produce muons that penetrate to the depth 
of the IceCube detector. Although few muons are produced in a gamma-ray shower rela- 
tive to a hadronic shower, a directional source of TeV photons can nevertheless produce a 
statistically significant excess over the large background of muons produced by cosmic ray 
primaries. To operate as a gamma ray observatory, IceCube records the directions and en- 
ergies of all "background" muons to create a muon sky map of the southern sky, something 
that no other experiment is capable of at this time. 

Clearly IceCube cannot match the instantaneous sensitivity of the new generation of 
atmospheric Qerenkov telescopes. The latter, however, are only capable of observing a 
several-degree patch of the sky on clear moonless nights. With a South Pole location, 
IceCube is unique in that it observes the same sky without interruption. It is also sensitive 
to sources overhead at the South Pole, a poorly studied portion of the southern sky. IceCube 
is the first large-scale detector in the Southern Hemisphere with the potential to detect TeV 
gamma-ray photons. The scientific potential of such all-sky TeV gamma-ray detectors has 
been clearly demonstrated by Milagro and Tibet AS7. 

Evaluating the capability of IceCube to observe Southern PeVatrons requires revisiting 
the estimates of the muon flux generated by photon showers. The sources of muons in a 
photon shower are threefold: the decay of mesons (pions and kaons) produced in photo- 
production events by shower photons on nuclei in the atmosphere, muon pair production 
in the electric field of a nucleus, and the production and decay of charm particles. At the 
energies considered in this paper muons of charm origin are few and their contribution is 
within the "errors" associated with the first two mechanisms [7] . Section [TTl discusses our 
analytical derivation of the muon flux created by a gamma-ray beam. Progress has been 
made possible by the appearance of simulation programs such as GEANT and CORSIKA. 
The EGS4 option in CORSIKA enables a full Monte Carlo simulation of the electromagnetic 
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component of showers and provides detailed information for all electromagnetic particles. 
Apart from standard photon interactions like Compton scattering and e^e~ pair production 
CORSIKA also includes direct ^"^^'-pair production and photonuclear reactions with nuclei 



in the atmosphere 



14l | . We can compare this with the earlier analytic shower calculations to 



converge on more reliable estimates. While the linear shower calculations match the shower 
calculations for the rate of muons of pionic origin, we found that this is not the case for 
muon pair production. The commonly followed procedure of substituting the electron by 
the muon mass in the expression for electron pair production results in underestimating 
the high-energy cross-section by a factor 2 and so our reevaluation of the muon pair cross 
section following the formula used by GEANT results in an enhanced production of high 
energy muons in photon showers compared to CORSIKA (version 6.900). 

Our model does suffer from certain systematic defects intrinsic to an analytical approach. 
For example, the flux we derive is necessarily time-averaged, resulting in our inability to 
predict the rate of muon bundles due to the occurrence of several muon production events 
within a single gamma shower. However, examination of CORSIKA event rates shows that 
the rate of multiple production events is small and does not affect the final significance of 
the signal. The strength of our approach lies in its transparency: all physical parameters 
and cross-sections can be controlled and easily changed if necessary. As a practical matter 
the code written for this work runs in < 10 s on a personal computer, much faster than a full 
Monte Carlo, making this approach the correct one for a comprehensive initial examination 
of the possibility of detecting gamma-induced showers and a good first approximation to 
the full simulation that will be required to perform an analysis of real experimental data. 

In Section |V] we discuss the primary sources of cosmic gamma-rays, in particular the 
PeVatrons producing pionic gamma rays, whose spectrum extends to several hundred TeV 
without cut-off, in interactions with the interstellar medium. By straightforward energetics 
arguments isl one anticipates the TeV flux from a source at a nominal distance of 1 kpc to 
be in the range 

Ti_iA in-i2 m-ii photons^ f W \ f n \ f d 
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E ™ (>1 TeV) = - 10'^" - 10 , 

dE AtkP \ cm^ s / \ 10^° erg/ \lcm~^/ \lkpc 

(1) 

Such sources must emerge in an all-sky TeV gamma ray survey performed with an instru- 
ment with the sensitivity of the Milagro experiment. Based on observations by Milagro and 
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H.E.S.S., it has been argued that one Pevatron, MGRO J1908+06, has hkely been identified 
among six candidates in the current sky map [2]. Our main conclusion (Figured]) is that a 
generic point-hke PeVatron with a normahzation at 1 TeV of 3 — 8 x 10^^^ TeV^^ cm^^ s~^ 
extending to 300 TeV may be observed with a significance of 3 — 5 a in 10 years by counting 
the total detected number of muons over the full energy range. We also emphasize the pos- 
sibility of exploiting the complementarity between all-sky and pointed telescopes that led to 
the final detailed observation and fiux measurements of J1908+06: once a significant excess 
has been observed in IceCube, a more sensitive pointed instrument can devote observation 
time to that specific region of the sky, with luck confirming the presence of a source and 
measuring the spectrum of gamma-rays. 
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FIG. 1: Differential flux as a function of spectral index from a point-like PeVatron required to create 
an excess over background with a statistical significance of 3 and 5 a (probabilities of 1.3 x 10~^ 
and 2.9 x 10"'', respectively) after 10 years. The dashed lines include an estimated contribution 
from kaon production and decay (Section IIC). 
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II. MUON FLUXES FROM GAMMA SHOWERS 



A. Gamma-ray cascades in the atmosphere 

The 7-rays that reach the Earth from the source will interact with an atom of the atmo- 
sphere producing an electron-positron pair, which will then via bremsstrahlung emit photons 
and create an electromagnetic cascade. Since the photons of the cascade can create muons 
that are detectable in IceCube, we must first find the spectrum of 7-rays at some slant depth 
t in the atmosphere 7 {E, t) created by the initial source flux. If our initial (differential in 
energy) spectrum can be represented by a power- law 7 (E, t = 0) = K^E-^'+'^\ cascade the- 



ory under Approximation A 16| tells us that the spectrum after cascading is a power-law of 
the same spectral index: 



A2 — Ai 



= 7(i?,t = 0)72(t), (2) 



(Jo + \l CTq + A2 

Ai^2 are the scale lengths of the shower growth and attenuation in the atmosphere and 
CTq ~ 7/9 is the probability of e^e~ pair production per radiation length. The A's are 



dependent on the integral spectral index s and are given in tabulated form in Ref. |16 |. 
Conventionally, t is in units of radiation (bremsstrahlung) length and so the parameters Aj 
and (To have units of (radiation length)"^. The most important special case is s = 1, because 
Ai(s = 1) = and A2(s = 1) < 0. As a result, an incident differential flux that goes as 
i?"^ will create a cascade that maximizes in the atmosphere after several radiation lengths, 
without any subsequent attenuation: 

7(E,t) ^ 0.5677 (£;,t = 0). (3) 

B. Pion production 

The principal channel for gamma-ray induced muons is by the photoproduction of pions 
on the nuclei of the atmosphere. Drees et al. derive a closed-form analytical solution for 
this flux using the linear cascade equations, assuming E""^ js]. In the following section, we 
generalize their method to other spectral indices, essentially starting with Eq. [2] rather than 
with Eq. [3] for the gamma-ray flux. Due to the explicit dependence on t in Eq. [21 the solution 



for a general spectral index cannot be expressed in a closed form and must be represented 
by an integral. However, this integral is not difficult to compute numerically and allows us 
to avoid the assumptions regarding depth that are required to have no final dependence on 
t, even for E""^. For a general derivation of the linear cascade equation and its simplification 
under Approximation A we recommend the discussion in Ref. [17 1. 

We begin with the ansatz that the pion spectrum (differential in energy) can be factorized 

as 

7r{E,t)=^{E,t = 0)7r2{E,t). (4) 

Where X-y(Tr)A and o"-y(^)A are respectively the interaction length and cross-section of pho- 
tons (pions) on nuclei of mass number A (in this case air nuclei with (A) = 14.5) and using 
the standard definitions of spectrum- weighted moment Zij, effective interaction length Aj, 
and pion decay length d^^it) [l7|, we can reduce the linear cascade equation 



to 



_XnA dj^it) _ ^ '' 
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The first term in Eq. [5] represents the loss of pions due to interactions and decay, and the 
second term represents the production of pions in photonuclear interactions. The production 
of lower-energy pions in interactions of pions with air, the third term in Eq. [5|, has the effect 
of reducing the loss of pions due to interactions and is taken into account by using the 
effective interaction length A^r in Eq. [H] which is larger than the normal interaction length 
Att. We assume that the inclusive differential cross-sections depend only on the fraction x 
of energy transferred from the parent particle (Feynman scaling), and that the interaction 
lengths A^(7r)A are constant over the energies considered here. 

We solve this differential equation in two regimes: high energy where pion interactions 
dominate over decay (A^ <^ d^^), and low energy where we can neglect pion interactions 
altogether (rf^ <^ A^). The energy scale is set by the pion decay energy constant = 115 
GeV in the decay length (i,r = EtcosO/tTr, where 6 is the zenith angle of the incident 
gamma-ray fiux. At high energy, where we can solve Eq. [6] by hand. 



At low energy, 



{cTQ + Al)((To + A2) 



A 



A2 — Ai 



gA2< g — t/A,r 



(ao + Ai)(Ai + ^) (cTo + A2)(A2 + 7^ 



(7) 



z^TT (o-Q + Ai)(ao + A2) 



'dt' (- 



.Alt' 



ctq + Ai ctq + A2 



(8) 



A7A A2 — Ai Jq \t 
where 6 = t/d-,,. 

For spectral indices such that Ai, A2 < the integral of Eq. M can be expressed as the 
lower incomplete gamma function. However, Ai > for spectra harder than E~'^, and 
moreover at low energies (5^1) and small depths the integrand will be so close to zero 
that the standard numerical implementations of the gamma function will fail. For these 
reasons we instead expand the exponential factor in the integrand as a series and perform 
the integration term-by-term: 



dt' 



t'V e^'*' 



100 



Ar' p 



(9) 



Due to the factorial in the denominator the series converges quickly and so we terminate 
the summation after 100 terms. The final pion flux is then an interpolation of the high- and 
low-energy limiting forms. However, since we do not have analytical expressions for both, 
smoothly transitioning from one limit to another is difficult and we take simply 



7i{E, t) = 7 {E, 0) min [vrf ^(t), ^^{E, t)] . (10) 

Comparing this to the smooth interpolation in Ref. jsl, we see that any discrepancy is at 
the transition from low- to high-energy behavior. This occurs at energies around e^r/ cos^ 
which is too low to produce detectable muons. Assuming no muon decay or energy loss in 
the atmosphere, standard 2-body decay kinematics 171] give the muon flux at the surface 



dN., 



dE 



dt B. 



E/r 



dE' n{E\t) 



(11) 



,0 JE (l-r)E' d^(t) 

where r = {m^/niT^Y and B^t^ = 1 is the number of muons per decaying pion. The upper 



cutoff on the depth integral tmax arises from the fact that the gamma-ray flux does not extend 
to arbitrarily high energies, so that at some depth in the cascade there should remain no 
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photons with sufficient energy to produce a muon of energy E. Following the 'Heitler model' 



argument of Ref. 



we take 

tmax — Ah m 



E ■ (12) 

where we define the effective Heitler cascade length = 8/7, the average of the 
bremsstrahlung and pair production lengths. For pionic decays the fraction of gamma-ray 
energy that goes into the final muon is (x)^^^ = 0.25. This form for tmax causes the final 
muon flux to depend approximately logarithmically on -Emax- For very low-energy muons, if 
we neglect muon decay the limiting factor is the physical extent of the atmosphere rather 
than the depletion of high-energy photons in the cascade, so tmax cannot be greater than 
(680/cos^^) g/cm^ (18.3/cos6' radiation lengths), the atmospheric depth of the South Pole 
ice surface. However, this depth is so large that it does not affect the muon flux that reaches 
the detector for any reasonable -Emax^ the muon energy that corresponds to tmax = 18.3 is 
given hj E = E^i^^ x (2.8 x 10^*^) and only the production of muons with less energy is 
suppressed due to the ice surface. Since the initial particle is a gamma-ray, taking Xh = 8/7 



underestimates the depth integral since the first electromagnetic cascade length dominates 
with Ag+e- = 9/7. Comparison with the Monte Carlo results of Ref. jo] and CORSIKA 
reveals that the analytic approach suffers from a systematic underestimation on the order 
of 20%, so in this work we multiply the muon fluxes from mesons by a factor 1.2. 

n 

We take the particle physics parameters used here from Ref. [8!]. We convert all lengths 
to radiation lengths, assuming the radiation length in air to be 37.1 g/cm^. We assume 
the photoproduction cross-section a^A = Aa^j^ = 14.5 x 100 /ib [5]. We obtain interaction 
lengths in air via 

A = 2.4 X 10^ (^i^^ [g/cm^] . (13) 

For the pion interactions A^^ = A^/(l — z^^t^) = 173g/cm^. Finally we assume that far from 
the resonance region equal numbers of positive, negative and neutral pions are produced 
which gives z^^^ = 2/3 for an incident E~'^ (s=l) spectrum, and that z^t^ and z,r7r vary 
sufficiently little that we can take the s = 1 value for all spectral indices of interest. 
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C. Kaon production 



At high energies, the production of muons from kaon decay becomes significant. A precise 
treatment of kaons would comphcate the determination of muon fluxes greatly, however. 
Kaons decay to pions, and pions can create kaons in interactions with atoms in the air, 
requiring not only an extra cascade equation for kaons but also additional source terms in 
the cascade equations that couple the meson fluxes. Fortunately, the coupling is weak: the 
z-factor for kaon production in pion interactions is of the order of 0.01 p^], and while the 
branching ratio for hadronic kaon decays is significant [iQ^, as is pion production in kaon 
interactions, the kaon flux is sufficiently small that its contribution to the pion flux can 
be neglected. As a result, we consider the fluxes of pions and kaons to be independent, 
neglecting the production both of pions by kaons and of kaons by pions. This allows us to 
re-use the pion formalism detailed above to estimate the kaon flux, substituting parameters 
where appropriate. We estimate the kaon production cross-section for energies far above 
threshold as a^N^K = cr-yN^ir {^tt/ttt^kY ■ This gives a^N^K ~ 8u6, a reasonable estimate 



since that is the order of the total cross-section near threshold |20l . I2l[ | and we expect a 
logarithmic increase in cross-section with energy to approximately counteract the fall in 
cross-section above the resonant threshold region. As with pions, we assume that 2/3 of 
the produced kaons are charged {z.^m^k = 2/3) and we neglect the neutral kaons. This 
results in underestimating the final muon flux since the neutral kaons typically decay either 
hadronically into pions or semileptonically into leptons and charged pions, but even were 
these to be treated here the uncertainties in the cross-section preclude estimating the kaon- 
induced muon flux with any great precision. We also neglect the hadronic decays of the 
charged kaon, considering only the 63.5% that decay to a muon and a muon neutrino: 
— s> /i^ + v^. Therefore B^k = 0.635. The scale energy of kaon decay in air is = 
850 GeV and we take the effective interaction length Aj^ = 200g/cm^. For kaon decays 
{x)^^fj, = 0.17 assuming an equal average fraction of gamma-ray energy into kaons as into 
pions. Due to the large uncertainties in the kaonic muon flux all event rates and statistical 
significances in this work are given both with and without the kaon contribution. 
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D. Muon pair production 



1. Cross-section 

One of the primary mechanisms by which an electromagnetic cascade proceeds is by the 
pair production of electrons by the shower photons. As a result, electron pair production has 



been understood since the earliest days of cosmic-ray physics [22f| but its generalization to 
more massive leptons has in many cases been flawed. In this section we describe the physics 
of pair production of massive leptons by an energetic photon, and discuss and correct the 
standard errors that have found their way into the literature and hence into standard software 
packages such as the Monte Carlo event generator CORSIKA (version 6.900). 

Pair production occurs when an incident photon interacts with a photon from the electric 
field of a nucleus, producing a pair of leptons. Since a single photon cannot produce a pair 
of massive particles without violating the conservation of 4-momentum, a second photon 
is necessary, transferring the required momentum from the nucleus. As the energy of the 
primary photon increases, the minimum momentum that the second photon must have 
for the pair production to occur decreases. As a result, the maximum distance from the 
nucleus at which the primary photon can still initiate pair production increases with energy 
and therefore we expect the cross-section to similarly increase. If the nucleus were a bare 
source of electric field we would expect the cross-section to increase without bound, as at 
extremely high energy the minimum momentum transfer necessary would be so small as 
to be possible at very large distances from the nucleus. However, nuclei are not electric 
monopoles out to arbitrarily large impact parameters. At a certain radius the nuclear 
electric field becomes 'screened' by the opposite electric field of the atomic electrons and at 
very large distances the atom is essentially electrically neutral. This decrease in effective 
electric charge beyond a certain distance from the nucleus cuts off the logarithmic increase 
of the cross-section and leads to the cross-section reaching an asymptotic maximum at high 
energy. Therefore, accurately predicting the cross-section requires estimating the maximum 
distance from the nucleus that pair production of a gamma-ray of energy into two leptons 
of energy and can occur, and then determining the effective nuclear charge visible at 
that distance and hence the probability of the emission of a virtual photon of the required 
momentum. The distance estimation is generally expressed by the screening parameter 
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S ~ ^z/-Rmax, essentially the ratio of the radius of the atomic cloud rz and the maximum 
impact parameter, i?max, to obtain the minimum necessary momentum transfer. From the 
conservation of energy and momentum and assuming energies large compared to the mass of 
the pair particles mi, we find the minimum necessary momentum transfer from the nucleus 
to be 

Therefore, the maximum distance from the nucleus at which this momentum transfer can 
occur is 

h 

-Rmax ~ • (15) 

Following the Thomas-Fermi model, the effective radius of an atom of atomic number Z is 



rz 



nip 



Z-'/\ (16) 



Therefore, the screening parameter is of the form 

S = "^^^^ (17) 

^i?„,ax 2E+E_(«meZV3)- ^ ^ 

Note that if in describing electron pair production we take m/ = rrie and cancel the mass 
factor in the denominator, the final expression will no longer be generalizable to other 
pair particle masses. 6 will be too small by a factor me/mi, resulting in the interaction 
reaching the full screening regime at overly low energies. For muon pairs on nitrogen, 6 = 1 
corresponds to a gamma-ray with energy ~ 100 TeV, which is therefore the approximate 
energy at which we expect screening to begin to saturate the cross-section. We see this by 
comparing in Fig. [2]curve 1, which uses the correct general expression for 6, to curve 2, which 
uses the reduced form of 6 correct for electron screening with the substitution m = m^. 

The second common error comes from the fact that the form factor of the atom must 
be integrated over the transferred momentum, where the upper limit is ~ mi. Therefore, 
the correct result for electron pair production is not generalizable to heavier particles and 
the limit at high energy will be too small by a factor of order ln(m//me) (curves 1 and 2 



in Fig. [2]) [71]. In this work, we take the muon pair production formula of Ref. j23|], adding 
also the inelastic contribution using the form factors from Ref. 2J| (curve 3 in Fig. [21). 
We note that the current version of the standard cosmic-ray Monte Carlo event generator 
CORSIKA (version 6.900) has been written to reproduce the cross-section of Ref. 25(] (curve 
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1), whose asymptotic value is a factor ~ 2 smaller than the correct value, and that this will 
beco„.ected.f„t.eve.i„„.y. 




FIG. 2: Muon pair production cross-section on air (Z=7.2): 1 has the correct screening but uses 



electron pair production atomic form factors [25| ; 2 is the Bethe-Heitler electron pair cross-section 
with m = ra^ ^\ ; 3 is the correct cross-section [23] with inelastic corrections 24 ] ; finally muon 



pair production via positron annihilation with atomic electrons (see Section [II El below) . 



2. Muon flux from pair production 

The spectrum of pair muons is derived from the cascade equation 
dn dxj{E^/x,t) 1 da^^f,^ f E 



. :^2^^^El^_:^llJpEEU^E^ = ^], (18) 
dt Jq X X^f, a^f, dx \ x J 

giving us, including a factor 2 for the two muons produced, and taking the cross-section in 
millibarns (mb) and depth t in radiation lengths. 
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dEf, 2.4 X 104 



Equation [19] assumes that both muons in the pair are separately measureable. In practice, 
however, detectors such as IceCube do not have the spatial resolution to distinguish between 
the two muons of the pair and will instead reconstruct a single muon with the total energy 
of the pair, the energy of the parent gamma-ray. The pair spectrum will therefore be 



dE 2.4 X 10^ Vg/cmV V mb 



dt-f2{t,s) (20) 



where 



dx—{x,E) (21) 



is the total cross-section for a gamma-ray of energy E to make a muon pair. 



E. Muon production via electrons 

Both processes outlined above, muon pair production and meson photoproduction, have 
analogues with an incident electron rather than an incident photon. It is possible for a 
cascade positron to annihilate an atomic electron, producing a muon pair + e~ —>■ fi~^ + 
fi~ and it is also possible for a cascade electron or positron to undergo an electronuclear 
interaction by producing a virtual photon, possibly resulting in the creation of a meson 

+ A ^ + 7* + A + v4 + tt/K. 

In practice, however, both these processes are negligible. The peak cross-section for muon 



pair production from electron pair annihilation is ~ 1 yub per electron [28| , giving ~ Z /ib 
per atom (Fig. [2]). Away from this peak, which occurs at a positron energy of 61 GeV, the 
cross-section falls off rapidly and is 0.2 Z fih at 500 GeV positron energy and is therefore 
essentially zero at energies that are high enough to produce muons with the > 600 GeV of 
energy that is needed to reach the detector through the ice overburden. By contrast, the 
cross-section for the electronuclear interaction is approximately 80% of the photonuclear 
interaction cross-section but the average fraction of energy that the virtual photon carries 
is only about 5% of the incident electron energy {29]. As a result, nuclear interactions 
of charged particles are an important energy loss mechanism without being a significant 
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source of secondary particles. A necessary caveat is that photonuclear interactions are not 
theoretically well-understood and the existing formulas for lepton-nucleus interactions are 
largely phenomenological descriptions designed to reproduce the data for muon energy loss in 
matter. The above numbers for electronuclear energy loss are derived from a muon energy 
loss formula using the substitution rrie, which is clearly a crude approximation. 

However, based on the example of muon pair production from a real photon, we do not 
expect our error due to this simple substitution to be greater than a factor ~ 2 and we can 
reasonably conclude that electronuclear interactions can be neglected in this analysis. 



III. MUON PROPAGATION AND ENERGY LOSS 

Neglecting muon energy loss and decay in the atmosphere, we treat muon energy loss in 
the ice above the detector by taking the standard average energy loss formula: 

' -a- hE (22) 



dx 



with a = 2.59 x lO^^TeV (g/cm^)-! and b = 3.63 x lO^^ {g/cm^)'^ in ice |30j. The general 
solution for final muon energy Ej at depth R given initial muon energy Ei is 

{Ef (R)) = {Ei + a/b) e-^^ - a/b (23) 
The flux at depth R is then, for Ei = {Ef + a/b)e^^ — a/b, 

{R, Ef) = — {Ei) e^^. (24) 



dEf' ' " dE, 

Assuming IceCube to be at a vertical depth of 1450 m of ice with density p = 0.92 g/cm^, 
we find that for a vertical muon to arrive at the detector with energy 100 GeV, the surface 
energy must be 607 GeV. Of course, as the range goes as d^ert/ cos 9, the necessary surface 
energy will increase with increasing zenith angle. 

In the case of muon pairs, where due to the lack of spatial resolution we consider the 
flux of pairs rather than the flux of muons, the energy loss of the pair will be the sum of 
the energy losses of the individual muons. Although the pair muons each take on average 
half of the gamma-ray energy, the distribution is fiat at low energies and has a minimum 
at a; = 0.5 at high energies. Therefore taking a simple average is a poor approximation of 
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the energy splitting and will overestimate the energy losses of the pair. We consider the 
average energy fraction taken by the higher-energy muon: at 1 TeV (a;)HE = 0.75 and varies 
by only 4% between 100 GeV and 1 PeV. The final energy of a muon pair of initial energy 
E is therefore taken to be the sum of the final energies of two muons with initial energies 
0.75 E and 0.25 

We assume that the angle between the directions of the muons and the direction of the 
parent gamma-ray is sufficiently small so that the distance between the muons on arrival 
at the detector is negligible compared to the spacing of the optical modules. While this is 
true at high energy, at low energies or large energy asymmetry between the muons one or 
both muons may miss the detector entirely. We neglect this possibility here as muons with 
a large angle to the gamma ray will be of low enough energy that they would not trigger 
the detector even if they were coUinear to the gamma-ray. 



IV. DETECTOR BACKGROUND AND EFFECTIVE AREA 



Finding sources of cosmic gamma-rays will depend on whether or not the signal in muons 
is sufficiently large to be detected against the background of cosmic-ray induced muons. 
As a result, we have to determine the spectrum and event rates of atmospheric muons in 
IceCube. We must also determine IceCube's response to downgoing high-energy muons to 
determine the signal event rates. To do this, we have used the atmospheric muon event rates 



in IceCube (3l| and have extracted effective areas from them. Given these effective areas 
we can then find the event rates from gamma-ray sources. We do not consider extensions 
to IceCube-80 such as the six-string cluster "DeepCore," as the additional optical modules 
are concentrated in the lower part of the detector and are therefore unlikely to significantly 
affect the effective area for low-energy downgoing muons. 

The integral event rate N{>E) is found by the convolution of the energy- dependent 
effective area with the muon flux at the detector, where T is the time interval: 

N,{>E) = tJ^ A,s{E,) ^ (E,) dE, (25) 
Therefore, given an atmospheric muon event rate we have to make an assumption about the 



cosmic-ray-induced muon flux to flnd t 



imation given in Chapter 6 of Ref. |l7| after propagation through the ice, and deconvolved 



le effective area. We have used the analytical approx- 
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the downgoing single muon effective area for three choices of zenith angle (Figure [3]). Due 
to the event rate being the output of a Monte Carlo simulation, it is not a smooth function, 
particularly at high energies where it is susceptible to fluctuations due to poor statistics and 
as a result, the derived effective area has corresponding unphysical "fluctuations." These do 
not affect our final conclusions as there are not sufficient events at those energies to affect the 
overall detectability of sources. At low energies the function matches our expectations for 
the downgoing muon effective area-an initial region of increasing area due to rising trigger 
efficiency followed by a region of approximately constant effective area, somewhat less than 
the standard 1 km^ due to the fact that IceCube's optical modules point down, reducing 
their sensitivity to light from downgoing muons. The 'trough' at intermediate energies is 
most likely due to the muon multiplicity of proton showers increasing with energy, which 
the time-averaged analytical approximation must neglect. 
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FIG. 3: Downgoing muon effective area for three zenith angles. Tfie fluctuations at high energy 
are due to fluctuations in the generated background event rate and do not affect the significance 
of signal event rates at lower energies. 

We assume that the muons contributing to the effective area have an angular resolution 
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of 0.5°, and find the radius of our circular search bin to determine the background flux 
magnitude from 



(bin radius) = 1.6 x a/ (angular resolution) ^ + (source radius)^. (26) 



This search bin will contain 72% of the total signal events 321]. The final significance of the 



excess above background in the search bin due to the signal is estimated using 

0.72 Nsi 

ria- = — (standard deviations). (27) 

V^bkg 

We use the number of background events A^bkg in the denominator rather than the more 
conservative A^bkg + ^sig, due to the fact that IceCube can measure the mean background 
in the bin at any time by averaging over the total events from the declination band of the 
sky at the source zenith. 



V. SOURCES OF COSMIC GAMMA-RAYS 



A. PeVatrons 



Point-like PeVatrons: A point source of cosmic rays with energies up to the knee 
at 3PeV is expected to be a source of gamma-rays with energies up to ~ E'kncc/lO = 
300 TeV. To model a generic PeVatron, we take a source with normalization = 
1 X 10^^^ TeV~^ cm~^ s^^ at 1 TeV and -Emax = 300 TeV. We assume an IceCube angu- 
lar resolution of 0.5° and a negligible source diameter. The muon flux from such a source, 
with incident gamma-ray spectral indices of -1.6, -2, and -2.6, is shown in Figs. HI [5l and El 
respectively. The event rates per year are shown in Fig. [71 The sources are assumed to be at 
a zenith angle of 30° but the final statistical significance (Fig. [8]) varies very little with zenith 
angles from 0° to 45°. This is due to the fact that while the mesonic muon fluxes (including 
the background) at high energy are proportional to 1/ cos^, the distance from the surface 
to the detector varies by the same factor and so the surface threshold is higher for sources 
with large zenith angle. Since the maximum significance is always at threshold. Fig. [1] shows 
the minimum normalization as a function of spectral index necessary to observe a point-like 
PeVatron at 5 a (P= 2.9 x 10^'') and 3a (P= 1.3 x 10^^) after 10 years' time without any 
cut on muon energy. From these numbers we determine that observing point-like PeVa- 
trons with spectra steeper than E^'^'^ will require normalizations at 1 TeV of approximately 
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3 — 8 X 10~^^ TeV~^ cm~^ (Figure[T]). Finally, in Figs. [9l[T0l andfTDwe show the fractional 
contribution to the total signal event rate from each component considered in Section HTl 

Extended PeVatrons: If a very strong extended source exists in the southern sky, such 
as a large molecular cloud complex powered by several cosmic ray beams, how large can 
it be to still be detectable by IceCube? We find that a circular source with normalization 
10~^ TeV^^ cm~^ s~^ at 1 TeV and spectral index E'"^ gives 5 a after 10 years with a diameter 
of 11.5° (not including kaons) or 13° (including kaons). 3o" is reached after 10 years by a 
source with diameter 19.5° (not including kaons) or 22° (including kaons). 



B. Vela Jr. 



The brightest TeV source in the Southern Hemisphere, Vela Jr. is a ~ 700-year-old 
supernova remnant almost directly in front of the Vela Pulsar, at a zenith angle of ~45°. 
H.E.S.S. observations of the remnant reveal a 2° diameter shell emitting gamma-rays with 
a normalization at 1 TeV of 1.89 x 10"^^ TeV^^ cm~^ s~^ and an spectral index 331]. 

H.E.S.S. does not measure above 20 TeV and since it is not known at what energy the 
spectrum cuts off, we show the differential flux in Fig. [12] with -Emax = 50 TeV and the 
significance after 10 years in Fig. [13] with -Emax = 50 TeV and 300 TeV. Although Vela Jr. is 
a strong source, its large diameter means that the background is large and that detection is 
unlikely. 



C. Other Sources 



We also considered the detection of extragalactic gamma-ray sources such as active galax- 
ies and gamma-ray bursts. Due to the showering, the flux of muons is strongly dependent 
on the upper energy cutoff of the gamma-ray spectrum. Unfortunately, the existence of 
extragalactic background photons means that an extragalactic source must be extremely 
close for photons of energy greater than ~ 10 TeV to reach the Earth 3J], so even sources 
that have large fluxes at TeV energies will not be detectable using shower-generated muons. 
We therefore do not consider extragalactic sources here. 

Finally, the possibility exists of detecting periodic gamma-ray sources such as LS5039, 
using the periodic modulation of the total event rate to establish the existence of a source. 
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FIG. 4: Differential flux of muons at the surface from a point-like PeVatron with E 




FIG. 5: Differential flux of muons at the surface from a point-like PeVatron with E' 
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FIG. 6: Differential flux of muons at the surface from a point-like PeVatron with E 




FIG. 7: Integral event rate per year from a point-like PeVatron for E E ^ and E 
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FIG. 8: Statistical significance after 10 years with an energy cut at Efj_ of excess muon events due 
to a point-like PeVatron with flux at 1 TeV of lO^^^ TeV^^ cm^^ g^i, for E~^-^, E'"^ and E''^-^, 
with kaons included (dotted) and without kaons (solid). 




FIG. 9: Fraction of total signal events from a PeVatron due to pionic, kaonic and pair muons for 
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FIG. 10: Fraction of total signal events from a PeVatron due to pionic, kaonic and pair muons for 
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FIG. 11: Fraction of total signal events from a PeVatron due to pionic, kaonic and pair muons for 
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FIG. 13: Statistical significance of Vela Jr. after 10 years for E'max = 50 TeV, 300 TeV, with kaons 
included (dotted) and without kaons (solid). 
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Unfortunately, the known sources are weak and have low enough cutoffs that there will likely 
be insufficient signal events per period to make a detection. It is also difficult to determine 
the significance of a periodic source without data. The signal is detected by searching for 
peaks in the Fourier transform of the total event rate, and the theoretical simulation of 
real data is done by adding random noise to the mean event rates. However, given the low 
signal event rates, it is easy for the noise to mask the source frequency, and changing the 
set of random numbers used greatly affects the height of power spectrum peaks. Hence, we 
cannot determine an "average" sensitivity to a hypothetical periodic source and must wait 
for actual data whose noise will either allow for the detection of a periodic source or will 
not. 

VI. CONCLUSIONS 

We have calculated the rates and sensitivities of IceCube as a TeV gamma ray observatory. 
IceCube, capable of detecting muons of energy ~0.1 TeV and above, can observe the presence 
of muons generated in multi-TeV gamma ray showers, and distinguish these events from 
background given a sufficiently bright source. While air Qerenkov telescopes are considerably 
more sensitive, they do not have the ability to monitor large portions of the sky continuously. 
Even though IceCube has a lower sensitivity to gamma-rays than Milagro, if a TeV-bright 
transient source occurs in the Southern Hemisphere, IceCube may be the only experiment 
capable of monitoring it. No additional hardware or software is needed in IceCube beyond 
its planned design. 

Our results indicate that a standard point source of very-high-energy photons is 
observable with IceCube provided its flux is greater than several 10~^^ (Figure [1]). By 
contrast, a large emission region composed of many individual PeVatrons is also observable 
given an extremely strong total flux, approximately 10~^ TeV~^ cm~^ s~^ at 1 TeV for a 
source of diameter greater than ~ 13°. 

In all cases, due to the large numbers of muons, the maximum statistical significance is 
obtained without any cut on muon energy. Since the flux of muons from pair production 
follows the parent gamma-ray spectrum rather than being a power steeper as in the case 
of mesonic muons, it might be expected that the greatest sensitivity would be at energies 
greater than several TeV due to the steepness of the atmospheric background. This does not 
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seem to be the case as the extremely large background of cosmic-ray muons is not sufficiently 
reduced at higher energies to compensate for the neglecting of many signal muons from pion 
decay. Therefore, the best search method seems to be to simply maximize the number of 
signal events by measuring the total number of down-going muons that trigger the detector. 
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